Boundary quantum critical phenomena with entanglement renormalization 
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We extend the formalism of entanglement renormalization to the study of boundary critical phe- 
nomena. The multi-scale entanglement renormalization ansatz (MERA), in its scale invariant ver- 
sion, offers a very compact approximation to quantum critical ground states. Here we show that, by 
adding a boundary to the scale invariant MERA, an accurate approximation to the critical ground 
state of an infinite chain with a boundary is obtained, from which one can extract boundary scal- 
ing operators and their scaling dimensions. Our construction, valid for arbitrary critical systems, 
produces an effective chain with explicit separation of energy scales that relates to Wilson's RG 
formulation of the Kondo problem. We test the approach by studying the quantum critical Ising 
model with free and fixed boundary conditions. 

PACS numbers: 03.67.-a, 05.50.+q, 11.25.Hf 



The multi-scale entanglement renormalization ansatz 
(MERA) 0, H is a tensor network introduced to effi- 
ciently represent ground states of quantum many-body 
systems on a lattice. It results naturally from a real 
space renormalization group (RG) transformation that 
employs unitary tensors [disentanglers) to remove short- 
range entanglement from the system, a process known 
as entanglement renormalization The use of disen- 
tanglers is a key difference between the MERA and the 
matrix product state (MPS) Q, another tensor network 
ansatz for quantum spin chains that is at the heart of 
Wilson's numerical RG (NRG) Q and of White's den- 
sity matrix RG (DMRG) [|. Both MPS and MERA can 
be used to describe ground states of translation invariant 
systems. In addition, however, and thanks to the dis- 
entanglers, the MERA can also incorporate scale invari- 
ance, thus becoming a natural ansatz to investigate fixed 
points of the RG flow. It has indeed been used to study 
non-critical RG fixed points, corresponding to systems 
with topological order [f|, and critical RG fixed points, 
corresponding to continuous quantum phase transitions 

The goal of this paper is to extend entanglement renor- 
malization to the study of boundary quantum critical 
phenomena. Given a semi-infinite ID lattice C at a (scale 
invariant) quantum critical point, we propose the use of 
the scale invariant MERA with a boundary as an ansatz 
for its ground state. The boundary of the MERA con- 
sists of a semi-infinite MPS that represents the ground 
state of an effective lattice £ obtained from C through 
an inhomogeneous coarse-graining, and where each site 
represents a different length scale. As in Wilson's for- 
mulation of the Kondo problem 4], lattice £ exhibits 
explicit scale separation of energies. This MPS is char- 
acterized by a single tensor from which boundary scaling 
operators and boundary scaling dimension, as well as the 
boundary contribution to the ground state energy and 
entanglement entropy can be extracted (see appendices) . 
We demonstrate our approach by studying a semi-infinite 
critical quantum Ising chain with free and fixed boundary 



conditions. 

Bulk MERA. — Recall that the scale invariant MERA 
for an infinite critical chain [J d, 0, H, B E3| is char- 
acterized by a unique pair of bulk tensors, namely a 
disentangler u and an isometry w, distributed in lay- 
ers according to Fig. QJi). A layer of disentanglers and 
isometries defines a real space RG transformation that 
can be used to coarse-grain the original lattice C, pro- 
ducing a sequence of increasingly coarse-grained lattices 
{£, £' , £" , ■ ■ ■}■ Under coarse-graining, a local operator 
o transforms according to the scaling superoperator >S, 

Fig.rjjii), 

o^o' ^o" ■■■ (1) 

A scaling operator 4> a is a special type of operator that, 
under coarse-graining, transforms into itself times some 
scaling factor. The scaling operators <j) a and scaling di- 
mensions A Q are obtained from the eigenvalue decompo- 
sition of the scaling superoperator S [g, [9j , 



S(4> a ) = X a 4>o 



A, 



log 3 A Q , 



(2) 



where the base 3 of the logarithm reflects the fact that the 
coarse-graining transformation maps three sites into one. 
Recall that from the scaling dimensions, which govern the 
decay of two-point correlators in the bulk, Fig. QJiii) , 



((/> a (r)(/) a (r')) 



\r — r'\ 



/|2A„ 



(3) 



one can extract the critical exponents of the model. More 
generally, the conformal data that characterizes the con- 
formal field theory (CFT) [ll[ associated to the critical 
chains can be extracted from the scale invariant MERA 
9] . [In this paper we only consider one-site scaling oper- 
ators, placed on special sites in the lattice, as opposed to 
more general two-site scaling operators, placed on arbi- 
trary sites (see Ref. |j[ for more details), since the latter 
are not essential to the present discussion on boundary 
critical phenomena.] 
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FIG. 1: (Color online) (i) The bulk (scale invariant) MERA is 
characterized by a bulk disentangler u and bulk isometry w. 
(ii) Scaling superoperator S for one-site local operators o, in 
terms of the bulk isometry w. (iii) The two-point correlator 
{4>a{r)(j> a (r')) for selected sites r and r can be computed by 
coarse-graining the lattice log 3 (|r — r'|) times, after which 
the two copies of 4>a are nearest neighbors and fuse into the 
identity with unit amplitude (by normalization of <f> a ). Each 
coarse-graning step multiplies <f> a by 3~ A °, resulting in the 
power-law decay of Eq. [3] (see Ref. Q for details) . 




FIG. 2: (Color online) (i) The boundary (scale invariant) 
MERA is characterized by the bulk tensors u and w and a 
boundary isometry m°. Copies of the boundary isometry w° 
are connected together, (ii) Scaling superoperator <S* in terms 
of the boundary isometry w°. (iii) The two-point correlator 
((^(0)0/3 (r)), for r = (3 t+1 - l)/2, is obtained by coarse- 
graining the lattice t « log 3 (r) times, after which and (f)f) 
are nearest neighbors and fuse into the identity with ampli- 
tude Co/3, see Eq. [5] 



Boundary MERA. — In order to represent the ground 
state of a semi-infinite, quantum critical spin chain, we 
propose to use the boundary MERA described in Fig. 
[DJi), which is made of semi-infinite layers of copies of the 
bulk disentangler u and bulk isometry w, together with 
copies of a boundary isometry to 4 placed at the ends of 
the layers. This ansatz defines a real space RG trans- 
formation that, away from the boundary, is identical to 
that of the bulk MERA, Eq. [TJ but such that a local 
operator o° at the boundary transforms according to the 
boundary scaling superoperator 5°, defined in terms of 
the boundary isometry w*, Fig. [2jii), 



(4) 



This allows us to identify a new set of scaling operators, 
namely the boundary scaling operators and scaling 
dimensions A*, which are obtained from the eigenvalue 
decomposition of <S°, 



(5) 



A correlator between a boundary scaling operator ^* and 
a bulk scaling operator (j>p, Fig. (Hiii), reads 



af3 



(6) 



r A»+A, 3 ■ 

A rather significant, particular case of the above is when 



we choose the boundary scaling operator to be the iden- 
tity operator I, which has vanishing scaling dimension, 



(Mr)) « % 



(7) 



Recall that in a critical system without a boundary, 
(<^s (»■■)} bum = f° r an y scann g operator (other than the 
identity). Eq. [7] tells us that in the boundary MERA, 
the presence of the boundary is felt everywhere in the 
bulk, with an intensity that decays as a power law with 
the distance to the boundary. This is precisely one of the 
trademarks of boundary critical systems, as described by 
a boundary CFT (BCFT) [UGll. 
Effective lattice. — Let 



H = h°(0) 



J2 h (r,r + 1) 



(8) 



be the Hamiltonian of the semi-infinite chain C. where h° 
is a boundary term and h(r, r+1) = h is a (constant) two- 
site interaction term corresponding to a critical fixed- 
point bulk Hamiltonian free of irrelevant operators and 
with vanishing ground state energy. That is, we assume 
that h is a (two-site 0) scaling operator with scaling 
dimension A = 2 (but see remark on transitional layers 
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FIG. 3: (Color online) (i) The effective lattice L is obtained 
by coarse-graining C in an inhomogeneous way. Hamilto- 
nian term K^' t+1 ^ involving 3 4 sites in L becomes the term 
3 1_t /i(t, t + 1) on C. (ii) Definition of h [1 ' 2] in terms of K [1 ' 2] , 
u and w. (iii) Applying one layer of coarse-graining to K^ 2 ' 3 ^ 
produces 3 _1 if ' 1 ' 2 ' . (iv) Copies of w are used to describe the 
ground state of H in £, Eq. 1101 (v) The ansatz is equivalent 
to a MPS. 



below), and thus normalized to (h) buyk = 0. Let 

Tt+X—l 

K [t,t+i]= J- h(r,r + l), rt = (3*-l)/2, (9) 

r=r t 

be an operator that collects the 3* two-site terms h in- 
cluded in the interval of sites [r t , ^t+i] of C. The bound- 
ary MERA defines an effective lattice C, with sites la- 
belled by t G {0, 1, 2, • ■ • oo}, that results from coarse- 
graining the original lattice C in an inhomogeneous way, 
such that site t G C corresponds to 0(3*) sites of C, see 
Fig. GD^i). Under the inhomogeneous coarse-graining, the 
original Hamiltonian becomes 

oo 

H = fc*(0,l) + j^A 1_ * h(t,t + l), (10) 
t=i 

where h°(Q, 1) corresponds to h°(0) + h(0, 1) in Eq. 
A is just the scaling factor (A = 3), and the two- 
site term A 1_t h(t, t + 1) results from coarse-graining 
j([t,t+i\_ p or instance, h(l, 2) comes from coarse-graining 
= h[l,2) + h(2,3) + /i(3,4) as in Fig. ^ii); the 

term 3 _1 /i(2,3) comes from coarse-graining if! 2 ' 3 ! into 
3-iif[i.2] j Fig jgj^ and then K [i,2] intQ ~ h ag be _ 

fore; more generally, K^' t+1 ^ is first coarse-grained into 




Distance from boundary, r Distance from boundary, r 



FIG. 4: (Color online) Left: numerical estimate (Z(r)) for free 
and fixed BC obtained with a boundary MERA. The exact 
solution approaches the bulk value 2/tt as « r _1 . Right: the 
error in (Z(r)} for free BC (very similar to that for fixed 
BC) shows that the non-vanishing expectation value of the 
bulk scaling operators, Eq. [7] is still accurately reproduced 
thousands of sites away from the boundary. 



3 i-t^[i,2] and tnen K [i,2] aga j n mto ^ To understand 
the origin of the scaling factor 3 _t , we notice that under 
coarse-graining in the bulk, K [*' t+1 ] (made of 3* terms h) 
becomes proportional to K t' -1 '*] (made of 3 t_1 terms h) 

A-[t.t+i] 3 -i K [t-i,t^ (t>l) (11) 

where the factor 3 _1 = 3~ 2 x 3 is due to the scaling 
dimension A = 2 of h (factor 3~ A = 3~ 2 ) and the fact 
that each term h in if I* -1 '*! comes from three terms h in 
K [t,t+i] (factor 3). Then, after t - 1 iterations of the RG 
transformation, K [*> t+1 ] indeed becomes 

K [tj+i] ^ &-t K [i,*\ ( t >i). (12) 

The effective Hamiltonian H is of the form derived by 
Wilson as part of his resolution of the Kondo problem 
4]. It describes a semi- infinite chain of sites that in- 
teract with nearest neighbors with the same interaction 
term h, which is multiplied by a factor A 1- ' decreasing 
exponentially fast with the distance to the boundary. No- 
tice, however, that while Wilson's derivation was for free 
fermions and the limit A — > 1 was eventually taken, here 
we started with a generic critical Hamiltonian and A re- 
mains fixed at A = 3. In a recent paper [l3|], Okunishi 
proposed an interesting generalization of Wilson's NRG 
approach that also considered a generic critical Hamil- 
tonian, although the scale A > 1 in Eq. [10] was simply 
introduced l by hand' as a regulator of the gapless spec- 
trum and made very close to 1 (e.g. A = 1.02) to avoid 
an 'undesired perturbation' in its degeneracy structure, 
and h was chosen to be simply equal to h. 

Optimization. — Given a critical Hamiltonian ii, Eq. 
|51 the bulk tensors u and w are computed using the op- 
timization algorithm for the bulk scale invariant MERA 
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FIG. 5: (Color online) A few boundary scaling dimensions, or- 
ganized in conformal towers [nj , for the quantum Ising model 
with free and fixed BC. We include bulk scaling dimensions 
for comparison. The primary fields identity I, spin a and 
energy e are identified with a circle. The boundary MERA 
accurately reproduces the smallest scaling dimensions of each 
conformal tower. 



discussed in Refs. [1, Then H is coarse-grained into 
the effective Hamiltonian H of Eq. [TUJ The boundary 
isometry w° is obtained with a simplified version (replac- 
ing the MERA with an MPS) of the energy minimization 
techniques used in Refs. 0, [l4[ for the scale invariant 
MERA. An important point is that, both for the bulk and 
boundary MERA, in practical simulations we consider a 
few (e.g. three in the example below) transitional layers 
made of tensors {u\, wi, Wi}, {u2, W2, w^}, etc. that are 
different from the fixed-point tensors {«,w,w*}. As dis- 
cussed in Ref. ({J [l4[ , these transitional layers allow us 
to modify the vector space dimension d of one site of the 
original model to some larger value x and, if (contrary to 
our assumption above) H contains irrelevant operators, 
significantly diminish their effect. 

Example. — In order to test the perf ormance of the 
approach, we use a boundary MERA [15| to approximate 
the ground state of the critical quantum Ising model on 
a semi-infinite chain, 

oo 

H lelns = V X(0) - ]T (X(r)X(r + 1) - Z(r)) (13) 

■r=0 

where X and Z are the Pauli matrices and the constant 
r\ determines whether the system has free BC (r\ = 0) 
or fixed BC (i] — ±1). Fig. 2] shows, through accurate 
estimates of the expectation value (Z(r)) for both free 
and fixed BC, that the boundary MERA indeed offers a 
good approximation to the ground state of the original 
lattice C, in spite of the fact that the ansatz uses, arbi- 
trarily close to the boundary, bulk tensors u and w that 
have been optimized in the absence of a boundary! Ta- 
ble U shows some boundary scaling dimensions obtained 
by diagonalizing the boundary scaling superoperator 5°. 
Finally, in Fig. [5l these scaling dimensions are organized 
according to the conformal towers of the primary fields 
predicted by BCFT 

In summary, we have extended the formalism of entan- 
glement renormalization to the study of a critical infinite 
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TABLE I: Some scaling dimensions for free and fixed BC. 



lattice with a boundary. The boundary (scale-invariant) 
MERA only depends on two tensors u and w that en- 
code bulk properties and a third tensor w° that encodes 
boundary properties. We have seen that this simple 
ansatz can accurately reproduce the expectation value 
of local observables near and away from the boundary, 
as well as the lower part of the spectrum of boundary 
scaling dimensions of the model. By considering a differ- 
ent boundary isometry at each level of coarse-graining, 
it is also possible to study boundary RG flows (e.g. from 
free BC to fixed BC). Our results provide a numerical 
route to the study of boundary conformal field theory [l2j 
that may find applications in several areas ranging from 
condensed matter physics (boundary critical behaviours 
and quantum impurity problems) to string theory (open 
strings and D-branes). 

Support from the Australian Research Council (APA, 
FF0668731, DP0878830) and the Spanish Ministerio de 
Ciencia e Innovation (RYC-2009-04318) is acknowledged. 



APPENDIX A: BOUNDARY-BULK FUSION 
AMPLITUDES 

In this appendix we explain how to extract the coeffi- 
cient Cap for the boundary-bulk correlator 

(ftWMr)) « ^gfe, (Al) 

and more generally the coefficient CZg that governs the 
fusion of a boundary scaling operator <f> a and a bulk scal- 
ing operator (jig sitting next to the boundary into another 
boundary scaling operator 

^(0)X^(1) ^ £^(0) (A2) 

7 

Recall that the boundary scaling operators <j) a are ob- 
tained from the eigenvalue decomposition of the bound- 
ary scaling superoperator 5°, Eq. [5J We can now con- 
sider the set of eigenoperators <p a of the superoperator 
5°* dual to the scaling superoperator <S°, 

= A*&, (A3) 

where — log 3 corresponds to the scaling dimension A* . 
By construction, the <^'s form a biorthonormal set with 

the 4> a 's, 

tr(O^) = 6 aB . (A4) 
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The boundary identity operator, 0q = F, has scaling di- 
mension Ao = and its counterpart 4>q is the one-site 
density matrix p° of the boundary. From p° we can com- 
pute for the two boundary sites, and then evaluate 

Cap as 



C Q/3 = tr[^ (2) (0,l)(C(0)® 



(A5) 



Fig. [5] shows this expression as well as the one corre- 
sponding to C%p, which involves (ftp, 0* and the 
boundary isometry w°. 
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FIG. 6: (i) Fusion amplitude C a p. (ii) Fusion amplitude CP* 
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FIG. 7: (Color online) (i) Density matrix p *- 2 - 1 in terms of p® 
and w° . (ii) Density matrix p in terms of p° and w° . The 
fixed-point density matrix p° has been denned in appendix A. 



For the quantum Ising model on a semi-infinite lat- 
tice we obtain, for free BC, a value E° = 0.I8I69023 
which is remarkably close to the exact solution -Ef xact = 
(1/2 - I/tt) = 0.I8I690I13... 0. For fixed boundary 
conditions we obtain a value E° — —0.45492968. Based 
upon the exact solution for finite chains of over a thou- 
sand sites, we estimate the error in this value is order 
1(T 6 . 



APPENDIX B: BOUNDARY GROUND STATE 
ENERGY 

In this appendix we explain how to determine the 
boundary contribution E° to the ground state energy, 
defined as the difference between the energy (H) of the 
semi-infinite chain and one half of the ground state en- 
ergy for an infinite chain, 

E» = (H)-- (Mf,r + l)U (Bl) 

r— — oo 

Since in this paper h is normalized such that it vanishes 
in the bulk (if this was not the case, we would simply 
choose h — (h) bulk as the new h), E® is the expectation 
value (H) of H in Eq. [5J In order to evaluate (H), we 
can coarse-grain the lattice C into the effective lattice £ 
and use that (H)c — (H)^ to write 

oo 

E? = (h*(p, l)) z + 3 1 - t (h(t,t+ . (B2) 

For any t > 1, (h(t 7 1 + 1))^ = tr(ph), where p is the two- 
site density matrix in C. Recall that, by construction, the 
semi-infinite MPS describing C is translation invariant for 
t = 1, 2, • • • , so that any two sites [t, t+l] (with t > 1) are 
indeed described by the same reduced density matrix p, 
which is obtained from the density matrix p° as indicated 
in Fig. [7] Since Ylt^i = 3/2, we arrive at 

E° = tv(p o{2) h°) + hi(ph), (B3) 
which is manifestly finite and can be easily computed. 



APPENDIX C: BOUNDARY ENTANGLEMENT 
ENTROPY 

In this section we discuss how to compute the bound- 
ary entropy from the boundary scale invariant MERA, 
from which we can extract Affleck-Ludwig's 1 around- state 
degeneracy 1 g associated to the boundary [17] . In the 
bulk, the von Neumann entropy of a block of L sites 
scales as [H, Q1, M, HI HI 

S bulk ( J L) = |log 2 (L)+c', (CI) 

where c is the central charge of the relevant CFT and 
c' is some non- universal constant, whereas for a block of 
L sites at the boundary of a semi-infinite chain, the von 
Neumann entropy scales as 0, [2l|, [22j 

S bmmAvT {L) = | log 2 (2L) + - + S b , (C2) 

where Sb = log 2 (g) is the boundary contribution to the 
entropy (l7| |. which depends on the boundary condition. 
As discussed in [23| , the entropy of a region of L = 3 T + 1 
sites in an infinite chain is the same as the entropy of a 
simpler, effective density matrix described in Fig. [5[i)- 
(ii). Similarly, the entropy of a block of L = (3 T + l)/2 
sites including the boundary of a semi-infinite chain is 
the same as the entropy of a simpler, effective density 
matrix described in Fig. [BJiii)-(iv). Fig. [5] shows the 
entropies S'bum and 5 bonnd „, for the Ising model. Noticing 
that St — S' boundary (L) — S' bulk (2L)/2, from which we can 
extract S b = 0.0007 w (that is, g = 1) for free BC and 
S b = -0.4992 w -1/2 (or g = 2" 1 / 2 ) for fixed BC. 
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FIG. 8: (i) Density matrix p of a block of L = 3 T + 1 sites in 
an infinite chain (without boundaries), for T = 2. The bulk 

(2) 

density matrix p|; ulk is the fixed-point, two-site density ma- 
trix of Ref. [9(.(ii) An effective density matrix with the same 
spectrum, obtained by removing disentanglers and isometries 
from p (which do not change the spectrum), (iii) Density 
matrix of a block of L = (3 T + l)/2 sites that includes the 
boundary of a semi-infinite chain, for T = 2. The fixed-point 
density matrix p° has been defined in appendix A. (iv) An 
effective density matrix with the same spectrum. 
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FIG. 9: Entanglement entropy scaling of bulk, and of free and 
fixed BC. Continuous lines represent the exact scaling. 




FIG. 10: (Color online) MERA for a finite system with open 
BC and critical bulk Hamiltonian, made of N — 4 x 3 T sites, 
T = 2. Most of this finite-size MERA is filled with copies 
of the bulk disentangler u and bulk isometry w optimized in 
an infinite system. The boundary isometries wf, and are 
optimized in a semi-infinite system. The top isometry wt is 
obtained by diagonalizing the coarse-grained Hamiltonian on 
the effective lattice Ct (£2 in this case). 



APPENDIX D: FINITE SYSTEM WITH OPEN 
BOUNDARY CONDITIONS 

In this appendix we discuss how to study the low en- 
ergy spectrum of a critical Hamiltonian on a finite system 
with two open boundaries. For this purpose, we consider 
a finite lattice C made of N sites with a Hamiltonian 

AT-2 

Hlr = hl(0) + Hr, r + 1) + h%{N - 1), (Dl) 

where h is a critical, fixed-point Hamiltonian term and 
the left and right open BC are specified by the boundary 
Hamiltonian terms h L and h^. 

For a system with = 4 x 3 T sites, we attempt to de- 
scribe the low energy subspace of Hlr with the MERA 
described in Fig. [TQ1 which consists of T layers of dis- 
entanglers and isometries. Each layer is filled with bulk 
tensors u and w, obtained by analysing an infinite system 
using the scale invariant MERA algorithm of Ref. Q . In 
adddition, each layer has a boundary isometry w L at its 
left end and a boundary isometry at its right end, 
obtained by analysing a semi-infinite system using the 
boundary scale invariant MERA described in this paper. 
Finally, the ansatz is completed with a top isometry wt, 
which has an open index labelling (an approximation to) 
different low energy eigenstates of Hlr- 

Thus, the low energy subsapce of Hlr is expressed in 
terms of five tensors {u, w, w L , w^, wt}- The bulk disen- 
tangler u and bulk isometry w encode the critical prop- 
erties of the bulk; the boundary isometries and 
encode the critical properties of each of the boundaries 
separately, and the top isometry wt is the only tensor 
that contains information about both boundary condi- 
tions, as well as about the finite size N of the system. 
This top isometry wt is simply determined by diago- 
nalizing the effective Hamiltonian H that results from 
coarse-graining the system with the T layers of disentan- 
glers and isometries. 

Fig. [TTJ shows the low energy spectrum of Hlr for 
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FIG. 11: Excitation spectra of the quantum Ising model on 
a finite lattice C of N = 4 x 3 5 = 972 sites (where each site 
is made of four spins, so that the system has 4 x 3 5 x 4 = 



spins) with open boundary conditions. The energy is 
expressed in units such that the gap between descendants is a 
multiple of unity. All inequivalent combinations of open BC 
are considered. The different open BC are (0) = free, (+) = 
fixed(up), ( — ) = fixed(down). 



all possible combinations of boundary conditions of the 
critical Ising model. In all cases, they agree with the 
predictions of BCFT 0. 



wd,l and wd,r, which are optimized by regarding them 
as a MPS- like structure representing the ground state of 
a coarse-grained lattice £., similarly as we did in the case 
of an open boundary, and which contain all the informa- 
tion about the defect, including defect scaling operators 
and their scaling dimensions. 

When tested with the quantum Ising model, we find 
[23| that the scale-invariant MERA with a fixed-point de- 
fect is capable of describing all conformal defects, namely 
the line of continuous Dirichlct (CD) conformal defects as 
well as the line of continuous Neumann (CN) conformal 
defects 0. 

Finally, we can use a slight modification of the above 
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APPENDIX E: DEFECTS AND INTERFACES 

Entanglement renormalization can also be used to 
study a defect in the bulk of a critical system on an infi- 
nite lattice C. 

For this situation we propose a scale invariant MERA 
with a special stripe of defect disentanglers and isome- 
tries, see Fig. HW i). For a fixed-point (i.e. scale invari- 
ant) defect, only three defect tensors are needed, namely 
a defect discntangler ud and left and right defect isomc- 
tries wd.l and wd,r- The scale-invariant MERA with a 
fixed-point defect is thus characterized by: (i) a pair of 
bulk tensors u and v, which are optimized in the absence 
of the defect and contain only information about the crit- 
ical theory in the bulk; and (ii) three defect tensors ud, 



FIG. 12: (i) Scale invariant MERA for a system with a defect, 
(ii) Scale invariant MERA to study the interface between two 
semi-infinite systems at criticality 



ansatz to connect two (different) critical semi-infinite lat- 
tices and study their interface. This is achieved by con- 
sidering two types of bulk tensors, namely left bulk ten- 
sors ul and wl for the system on the left and right bulk 
tensors ur and wr for the system on the right, in addition 
to interface tensors itj, wi y L and wj,_r, see Fig. [121 ii). For 
instance, this ansatz can be used to study fixed-point in- 
terfaces between a semi-infinite quantum Ising chain and 
a semi-infinite quantum XX chain [23| . 
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